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Abstract. Some features of the global entanglement of a composed quantum system can be 
quantified in terms of the purity of a balanced bipartition, made up of half of its subsystems. 
For the given bipartition, purity can always be minimized by taking a suitable (pure) state. 
When many bipartitions are considered, the requirement that purity be minimal for all 
bipartitions can engender conflicts and frustration arises. This unearths an interesting link 
between frustration and multipartite entanglement, defined as the average purity over all 
(balanced) bipartitions. 
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1. Introduction 

Frustration in humans and animals arises from unfulfilled needs. Freud related frustration 
to goal attainment and identified inhibiting conditions that hinder the realization of a given 
objective flU. In the psychological literature one can find many diverse definitions, but 
roughly speaking, a situation is defined as frustrating when a physical, social, conceptual 
or environmental obstacle prevents the satisfaction of a desire 0. Interestingly, definitions 
of frustration have appeared even in the jurisdictional literature and appear to be related to an 
increased incidence of parties seeking to be excused from performance of their contractual 
obligations 0. There, "Frustration occurs whenever the law recognises that without default 
of either party a contractual obligation has become incapable of being performed because the 
circumstances in which performance is called for would render it a thing radically different 
from that which was undertaken by the contract.... It was not this I promised to do." JH 

In physics, this concept must be mathematized. A paradigmatic example flU is that of 
three characters, A, B and C, who are not good friends and do not want to share a room. 
However, there are only two available rooms, so that at least two of them will have to stay 
together. Their needs will therefore not be fulfilled and frustration will arise. A schematic 
mathematical description of this phenomenon consists in assigning a "coupling" constant 
to each couple (i, k), with i,k = A, B, C: = +1(— 1) if i and k (do not) like to share 
a room. Each character is then assigned a dichotomic variable Si = +1(— 1) if i is in the 
first (second) room. The key ingredient is the definition of a cost function that quantifies the 
amount of "discomfort" (unfulfilled needs) of our three characters. This can be easily done: 



Each addendum in the summation can take only two values, ±1. However, although each 
single addendum can be made equal to —1 (separate rooms, minimum cost and no discomfort 
for the given couple), their sum, in the best case, is if m in = — 1, which is larger than the sum 
of the three minima, —3. At least two characters will have to share a room and frustration 
arises. The situation becomes more complicated (and interesting) when more characters are 
involved and the coupling constants in (OQ) are, e.g., statistically distributed. 

The above description of frustration, in terms of a cost function, applies to a classical 
physical system. Interestingly, there is a frustration associated with quantum entanglement in 
many body systems. The study of this problem will be the object of the present investigation. 

Entanglement is a very characteristic trait of quantum mechanics, that was identified at 
the dawn of the theory Ed HI, has no analogue in classical physics BH and came recently 
to be viewed as a resource in quantum information science ifTOl [TT| . When the system is 
bipartite, its entanglement can be unambiguously quantified in terms of the von Neumann 
entropy or the entanglement of formation lfT2l [T3ll . Difficulties arise, however, when one 
endeavours to define multipartite entanglement [[141 CGI EOS IH [EH . The main roadblock is 




(1) 



The goal is to minimize this cost. In our case 



1, Vz, k, so that 



H — SaSb + SaSc + SbSc- 



(2) 
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due to the fact that states endowed with large entanglement typically involve exponentially 
many coefficients and cannot be quantified with a few measures. This interesting feature of 
multipartite entanglement, already alluded to in Ref. |[T9ll . motivated us to look for a statistical 
approach ll20ll . based on a characterization of entanglement that makes use of the probability 
density function of the entanglement of a subsystem over all (balanced) bipartitions of the total 
system ||2"T|| . A state has a large multipartite entanglement if its average bipartite entanglement 
is large (and possibly also largely independent of the bipartition). 

Maximally multipartite entangled states (MMES) ||22ll23l are states whose entanglement 
is maximal for every (balanced) bipartition. The study of MMES has brought to light the 
presence of frustration in the system, highlighting the complexity inherent in the phenomenon 
of multipartite entanglement. Frustration in MMES is due to a "competition" among 
biparititions and the impossibility of fulfilling the requirement of maximal entanglement for 
all of them, given the quantum state |[22l[20l . The links between entanglement and frustration 
were also investigated in Refs. Il24ll25ll26l . 

This paper is organized as follows. We introduce notation and define maximally bipartite 
and maximally multipartite entangled states in Sec.|2] We numerically investigate these states 
and show that multipartite entanglement is a complex phenomenon and exhibits frustration, 
whose features are studied in Sec. [3] Section @] contains our conclusions and an outlook. 

2. From bipartite to multipartite entanglement 

The notion of MMES was originally introduced for qubits [|22l and then extended to 
continuous variable systems ll27l l28l . Here we follow ll28ll and give a system-independent 
formulation. Consider a system composed of n identical (but distinguishable) subsystems. Its 
Hilbert space H = Hs, with Hs '■= ®ies ^ aiK * S = {!> 2, • • • , n}, is the tensor product of 
the Hilbert spaces of its elementary constituents hj ~ h. Examples range from qubits, where 
f) = C 2 , to continuous variables systems, where f) = L 2 (R). We will denote a bipartition of 
system S by the pair (A, A), where A C S, A = S\A and 1 < ua < ua, with ha = \A\, the 
cardinality of party A (n A + ua = n). At the level of Hilbert spaces we get 



Let the total system be in a pure state \ip) G H, which is the only case we will consider 
henceforth. The amount of entanglement between party A and party A can be quantified, for 
instance, in terms of the purity 



H = H A ® U a . 



(3) 



7T A = tr(p 2 A ) 
of the reduced density matrix of party A, 

PA = tT HA M)(ll>\). 

Purity ranges between 



(4) 



(5) 




(6) 



where 



mm 



(dim Ha) 



(dim h) 



(7) 
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with the stipulation that l/oo = 0. The upper bound 1 is attained by unentangled, factorized 
states = \4>)a ® \x)a (according to the given bipartition). When dimt) < oo, the lower 
bound, that depends only on the number of elements tia composing party A, is attained by 
maximally bipartite entangled states, whose reduced density matrix is a completely mixed 
state 

Pa = *&1ha (8) 
where 7r™f n is defined in Eq.© and 1- Ha is the identity operator on Ha- This property is valid 
at fixed bipartition (A, A); we now try and extend it to more bipartitions. 

Consider the average purity ("potential of multipartite entanglement") ||29ll22l 

*8Lm)=nu = (l ) E ( 9 ) 

V A J \A\=n A 

where E denotes the expectation value, the combinatorial coefficient is the number of 
bipartitions, \A\ is the cardinality of A and the sum is over balanced bipartitions tia = 
[n/2], where [•] denotes the integer part. The quantity 7Tme measures the average bipartite 
entanglement over all possible balanced bipartitions and inherits the bounds © 

< *8Lm < i. do) 

A maximally multipartite entangled state (MMES) [|22| | ip) is a minimizer of 7Tme, 

«£L(\<p)) = 4 n \ (id 

with = min{7rS(|^)) | |^> G H s , = 1}- 

The meaning of this definition is clear: most measures of bipartite entanglement (for pure 
states) exploit the fact that when a pure quantum state is entangled, its constituents are in a 
mixed state. We are simply generalizing the above distinctive trait to the case of multipartite 
entanglement, by requiring that this feature be valid for all bipartitions. The density matrix of 
each subsystem A C S of a MMES is as mixed as possible (given the constraint that the total 
system is in a pure state), so that the information contained in a MMES is as distributed as 
possible. The average purity introduced in Eq. © is related to the average linear entropy ll29l 
and extends ideas put forward in |[T7ll30l . 

We shall say that a MMES is perfect when the lower bound (ITOT) is saturated 

Et ] =ram{^l} = ^. (12) 

It is immediate to see that a necessary and sufficient condition for a state to be a MMES is 
to be maximally entangled with respect to balanced bipartitions, i.e. those with ua = [n/2]. 
Since this is a very strong requirement, perfect MMES may not exist for n > 2 (when n = 2 
the above equation can be trivially satisfied) and the set of perfect MMES can be empty. 

In the best of all possible worlds one can still seek for the (nonempty) class of states that 
better approximate perfect MMESs, that is states with minimal average purity. We shall say 
that a MMES is uniformly optimal when its distribution of entanglement is as fair as possible, 
namely when the variance vanishes: 

1/2 



a 



W = E 



(n) 

71 A - VT ME 



0, (13) 
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where the expectation is taken according to the same distribution as in Eq. © (all balanced 
bipartitions). Of course, a perfect MMES is optimal. It is not obvious that uniformly optimal 
non-perfect MMES exist. 

The very fact that perfect MMES may not exist is a symptom of frustration. We 
emphasize that this frustration is a consequence of the conflicting requirements that 
entanglement be maximal for all possible bipartitions of the system. 

2.1. Qubits and the symptoms of frustration 

For qubits the total Hilbert space is Hs = (C 2 )® n and factorizes into Hs = H A <8> H A , 
with H A = (C 2 ) 0nA , of dimensions N A = 2 nA and N A = 2 n *, respectively (N A N A = N). 
Equations © and rfTQfr -dTTb read 



respectively. 

For small values of n one can tackle the minimization problem (fTTI) both analytically and 
numerically. For n = 2, 3, 5, 6 the average purity saturates its minimum in (fl~5l) : this means 
that purity is minimal for all balanced bipartitions. In this case the MMES is perfect. 

For n = 2 (perfect) MMES are Bell states up to local unitary transformations, while 
for n = 3 they are equivalent to the GHZ states BTI . For n = 4 one numerically obtains 
E [ q ] = minvr^E = 1/3 > 1/4 = 1/N A [|22l [32l [33l |3l- For n = 5 and 6 one can find 
several examples of perfect MESS (22l[2l. The case n = 7 is still open, our best estimate 
being eP ~ 0.13387 > 1/8 = 1/N A . Most interestingly, perfect MMES do not exist for 
n > 8 11291 . These findings are summarized in Table Q] (left column) and bring to light the 
intriguing feature of multipartite entanglement we are interested in: since the minimum 1/N A 
in Eq. (fl~5l) cannot be saturated, the value of n A must be larger for some bipartitions A. We 
view this "competition" among different bipartitions as a phenomenon of frustration: it is 
already present for n as small as 4. This frustration is the main reason for the difficulties 
one encounters in minimizing 7Tme in ©. Notice that the dimension of Hs is N = 2 n and 
the number of partitions scales like N. We therefore need to define a viable strategy for the 
characterization of the frustration in MMES, even for relatively small values of n. 

2.2. Continuous variables and further symptoms of frustration 

For continuous variables we have Hs = (L 2 (M))' x,n with dim L 2 (M) = oo. As a consequence, 
the lower bound 7r™f n = in Eq. © is not attained by any state. Therefore, strictly 
speaking, in this situation there do not even exist maximally bipartite entangled states, but 
only states that approximate them. This inconvenience can be overcome by introducing 
physical constraints related to the limited amount of resources that one has in real life. This 
reduces the set of possible states and induces one to reformulate the question in the form: 
what are the physical minimizers of ©, namely the states that minimize © and belong to the 



l/N A <71 A <1, 

i/n a < < K^im < i 



N A = 2 [n/2 \ 



(14) 
(15) 
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Table 1. Comparison between qubit and Gaussian maximally multipartite entangled states for 
different number n of subsystems. 



n 


qubit perfect MMES 


Gaussian perfect MMES 


2,3 


yes 


yes 


4 


no 


no 


5,6 


yes 


no, but uniformly optimal* 


7 


no* 


no 


> 8 


no 


no 



*numerical evidence 



set C of physically constrained states? In sensible situations, e.g. when one considers states 
with bounded energy and bounded number of particles, the purity lower bound 

<t C = infW, m eC}> <fc (16) 

is no longer zero and is attained by a class of minimizers, namely the maximally bipartite 
entangled states. If this is the case, we can also consider multipartite entanglement and ask 
whether there exist states in C that are maximally entangled for every bipartition (A, A), and 
therefore satisfy the extremal property 

*a = KSf (17) 
for every subsystem A C S with ha = \A\ < n/2. In analogy with the discrete variable 
situation, where dim f) < oo and C = Ti, we will call a state that satisfies (flTT) a perfect 
MMES (subordinate to the constraint C). 

Since, once again, the requirement (fTTI) is very strong, the answer to this quest can be 
negative for n > 2 (again, when n = 2 it is trivially satisfied) and the set of perfect MMES can 
be empty. In the best of the best of all possible worlds one can still seek for the (nonempty) 
class of states that better approximate perfect MMESs, that is states with minimal average 
purity. In conclusion, by definition a MMES is a state that belongs to C and minimizes the 
potential of multipartite entanglement ©. Obviously, when 

P W,C _ •„ r(n) x _ [n/2],C m 

E - mm{vr ME ) - vr min (18) 
there is no frustration and the MMESs are perfect. Eventually, we will consider the limit 

Let us consider the quantum state \il>( n )) of n identical bosonic oscillators with 
(adimensional) canonical variables {qk,Pk}k=i,...n an d un it frequency (set h = 1). An 
analogous description of the system can be given in terms of the Wigner function on the 
n-mode phase space 



W {n) (q,p) = J d n y(q-y\tp {n) )(tlj {n) \q + y)e 2 ™y p , (19) 

where q = (q 1} ... q n ), p = ( Pl , . . . p n ), y = (y u . . . y n ) £ W 1 , and we have denoted by 

\q±y) = ® n k=l \q k ±y k ) (20) 
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the generalized position eigenstates. By definition, Gaussian states 11351 [361 are those 
described by a Gaussian Wigner function. Introducing the phase-space coordinate vector 
X = (Xl, . . . X 2n ) = {qi,Pi, ■ ■ ■ Qn,Pn), a Gaussian state has a Wigner function of the 
following form: 

W in) (X) 



(27r)" v /det(V) 
1 



x exp 



2 {X - x )v-\x - X ) r 



(21) 



where X = (X) = J XW( n )(X)d 2n X, is the vector of first moments, and V is the 2n x 2n 
covariance matrix, whose elements are 

Y lm = ((X, - (X»(X m - (X m ))). (22) 

For Gaussian states, purity is a function of the "sub"determinant of the covariance matrix 

tta = ; = , (23) 

2^v^t(Vl)' 

Va being the square submatrix defined by the indices pertaining to bypartition A. Clearly, 

7T A < 1. 

The results of the search for perfect and uniformly optimal MMES with Gaussian states 
are summarized in the right column of Table Q] ll28l . There are curious analogies and 
differences with qubit MMES. In particular, for n > 8 perfect MMES do not exist in both 
scenarios. Actually, for Gaussian states, frustration is present for n > 4; for the "special" 
integers n = 5, 6 we notice that both for two-level and continuous variables systems the 
variance of the distribution of entanglement goes to zero; on the other hand, in the former 
case one can find perfect MMES, in the latter case MMES are uniformly optimal but not 
perfect. 



3. Scrutinizing Frustration 

We now turn to the detailed study of the structure of frustration. For the sake of concreteness, 
we shall first focus on qubits. Let us start from a few preliminary remarks. We observed that 
it is always possible to saturate the lower bound in (fT4l) 

7T A = l/N A (24) 

for a given balanced bipartition (A, A) with Na = 2' n / 2 l However, in order to saturate the 
lower bound in (fl~5l) 

/4 n) = 1/N A , (25) 

condition (|24|) must be valid for every bipartition in the average ©. As we mentioned in Sec. 
12. 1[ this requirement can be satisfied only for very few "special" values of n (n — 2, 3, 5 and 
6, see Table [I])- For all other values of n this is impossible: different bipartitions "compete" 
with each other, and the minimum of is strictly larger than 1 /Na- 
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Figure 1. (Color online) Qubits: average purity \i of the typical states (squares), average purity 
7Tme(SDC) computed according to Scott's procedure |29l with extremal additive self-dual 
codes (full circles), our best numerical estimate Eq (num) for the minimum of 7Tme (diamonds) 
and lower bound LB = 1/Na (triangles) vs the number of spins n. Note that for n = 7 the 
full circle and the diamond do not coincide. The scale on the ordinates is logarithmic. 



It is interesting to look at this phenomenon in more detail. Let us recall that for typical 
states []37l|38l[39l|40l|4Tl|2llll2l the distribution of purity over balanced bipartitions has mean 



Figure CD displays the average purity /i of typical states [Eq. (|26|)1. the average purity of 
extremal additive self-dual codes states, computed according to Scott's procedure If29l . our 
best numerical estimate for the minimum of Eq, and the lower bound 1/Na- All these 
quantities exponentially vanish as n — > oo. Scott's states give an upper bound for the minimal 
average purity when n > 8, where the numerical simulations become very time consuming. 
In particular, we notice that for 3 < n < 6 and n = 8 the numerical values of Eq coincide with 
the results obtained using extremal additive self-dual codes ||29ll . For n = 7 the optimization 
algorithm reaches a lower value. For n > 8 our numerical data do not enable us to draw any 
conclusions. 

Figure [2] displays the (normalized) difference between the minimum average purity, 
computed according to Scott's extremal additive self-dual codes and/or our best numerical 
estimate, and the lower bound 1/Na m Eq. (T25T) . This difference is an upper bound to the 
frustration ratio 

F (n) = E n) - 1 / N A Na = 2 M (27) 

^0 

that can be viewed as the "amount of frustration" in the system. We notice the very different 
behavior between odd and even values of n. In the former case the amount of frustration 
increases with the size of the system. On the contrary, in the latter case, the behavior is 
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not monotonic. It would be of great interest to understand how this quantity behaves in the 
thermodynamical limit, but our data do not enable us to draw any clear-cut conclusions. 

Let us extend these considerations to the continuous variables scenario. In order to 
measure the amount of frustration (for states belonging to the constrained set C) we define 
a more general frustration ratio 

tf(n),C _ [n/2],C 

F {n),C _ ffO ^min (2 Q\ 

^ ~ F (n),C 

and eventually take the limit C — > H, where both the numerator denominator can vanish. 
[Notice that (1271) is a specialization of the quantity in d28l to the qubit case, i.e. C = H = 
(C 2 )®".] As a constraint we fix the value J\f of the average number of excitations per mode, 
namely 

C = {ip E H, ip Gaussian, ^ + Pk ^ < M + \, 1 < k < n}. (29) 

_ 



The ideal lower bound is given by 11281 

tP!W = - (30) 

and represents the purity of a Gaussian thermal state. Incidentally, we notice that, for 
Af = 1/2, Eq. (l30l) reproduces the lower bound of Eq. (TT5T) i. e. the case of qubits. On 
the other hand, J\f = corresponds to the case of completely separable states. In Fig. [3] we 
plot the frustration ratio (l28l as a function of the number of modes. Each point has been 
numerically obtained by relaxing the energy constraint (J\f — > +oo) until the ratio F^' c has 
reached a saturation value ll28l . This corresponds to the limit C — > H. 



3.1. The structure of frustration 

In order to try and understand the underlying structure of this frustration, we focus again on 
qubits and analyze the behavior of the minimum value of the average purity when one requires 
the condition ((24)) for an increasing number of bipartitions. We proceed as follows: we set 
n = 7, so that the total number of balanced bipartitions is K = 35, and add bipartitions one 
by one, by choosing them so that condition (124)) be valid, as far as this is possible. When 
(1241 becomes impossible to satisfy, we require that the average purity be minimal. We plot 
the minimum average purity as a function of k and a = k/K in Fig. |4] . One observes that 
it is possible to saturate the minimum 1/Na = 1/8 = 0.125 up to k = 32 well chosen 
partitions. For k = 33 all bipartitions yield the minimum 1/8, except the last one, that yields 
1/4: frustration appears. For k = 34, two bipartitions yield 1/4. For k = K = 35 purity is 
larger than 1/8 for all bipartitions. Notice that the solution with 32 bipartitions at 1/8 and the 
remaining three at 1/4 corresponds to Scott's extremal additive self-dual code and would yield 
a higher average. The distribution of purity for an increasing number of partitions is shown in 
Fig. El 

The case n = 8 is slightly different: see Fig. |6l In this case K = 35 again and there 
is no frustration up to k = 28 bipartitions (a = 0.8), if properly chosen (all of them with a 




Figure 2. (Color online) Qubits: frustration ratio: (normalized) distance between the average 
purity 7Tme(SDC), according to Scott's extremal additive self-dual codes (full circles), or our 
best numerical estimate for the minimum Eq (num.) (squares), and the lower bound LB=1 /Na, 
vs the number of spins. Notice the difference between even and odd n. 




4 5 6 7 8 9 

Figure 3. (Color online) Gaussian states: frustration ratio for and odd (rectangle) and even 
(circle) number of modes. As explained in the text, each point has been obtained in the 
saturation regime C — > TL. 



purity 1/16). When k is further increased, it is no longer possible to reach the lower bound: 
for 29 < k < 32 the new bipartitions have purity 1/8. Finally, for 33 < k < 35 = K, the 
new bipartitions have purity 1/4. 

Another useful test is the extraction of k randomly selected bipartitions and the 
successive minimization of the average purity. This is a typical test in frustrated systems, 
e.g. in random ll43l and Bethe lattices ll44l . In this way one checks the onset of frustration 
independently of the particular choice of the sequence of bipartitions. In Fig. [7] we plot the 
dependence on k and a. Each point corresponds to the extraction of a number of k bipartitions 
ranging from few tens to a few hundreds. Frustration appears at rather large values of k, 
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qualitatively confirming the result shown in Fig. @J Moreover we notice that for smaller k 
it is sometimes difficult to reach the minimum. This could be an indicator that for a small 
number of bipartitions there is a large number of local minima in the energy landscape, that 
"traps" the numerical procedure. Notice that the curve in Fig.fTJshould monotonically increase 
as a function of k, so that all deviations from monotonicity are ascribable to the numerical 
procedure, and are a consequence of the fact that for different values of k, in each run of the 
simulation, the subset of extracted bipartitions is uncorrelated to the set used for the preceding 
values of k. 

Although the results of this section are not conclusive, they provide a clear picture of the 
relationship between entanglement and frustration. The latter tends to grow with the size of 
the system (Figs. [2] and [3]) and it is difficult to study, at least for small values of n, because 
it suddenly appears at the last few bipartitions (Fig. |4]). One estimates, from the results for 
n = 7 qubits in Fig. [71 an average fraction of frustrated bipartitions 1 — a c ~ 34%, a c being 
a critical ratio. 

A posteriori, it is not surprising that multipartite entanglement, being a complex 
phenomenon, exhibits frustration. It would be of great interest to understand what happens 
for larger values of n. Different scenarios are possible, according to the mutual interplay 
between the quantities F and a c . In particular, the amount of frustration F shows a tendency 
to increase with n, for both qubit and Gaussian states. This could be ascribable to a decrease 
of a c , corresponding to an increasing fraction of frustrated partitions, or to a constant (or even 
increasing) a c , corresponding to a constant (or decreasing) fraction of increasingly frustrated 
bipartitions. 

4. Concluding remarks 

One important property that we have not investigated here and that is often used to 
characterize multipartite entanglement is the so-called monogamy of entanglement lfi~4l 
l45l . that essentially states that entanglement cannot be freely shared among the parties. 
Interestingly, although monogamy is a typical property of multipartite entanglement, it is 
expressed in terms of a bound on a sum of bipartite entanglement measures. This is 
reminiscent of the approach taken in this paper. The curious fact that bipartite sharing of 
entanglement is bounded might have interesting consequences in the present context. It would 
be worth understanding whether monogamy of entanglement generates frustration. 

Two crucial issues must be elucidated. First, the striking similarities and small 
differences between qubits and Gaussian MMES: see Table [Hand compare Figs. [2] and |3j 
Second, the features of MMES for n — > oo. Finally, we think that the characterization 
of multipartite entanglement investigated here can be important for the analysis of the 
entanglement features of many-body systems, such as spin systems and systems close to 
criticality. 
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Figure 4. (Color online) Minimum average purity for n = 7 and an increasing number k of 
bipartitions. Until k = 32 no frustration appears. For k = 33 one partition has purity tt = 1/4. 
For k = 34, two bipartitions yield tt = 1/4. For k = K = 35, it > 1/8 for all bipartitions. 
a = k/K. See Fig. 
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Figure 5. (Color online) Distribution of purity for n = 7 and an increasing number k of 
bipartitions. See Fig. [4] 
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Figure 6. (Color online) Minimum average purity for n — 8 and an increasing number k of 
bipartitions. Until k = 28 no frustration appears, a = k/K. 



k 

5 10 15 20 25 30 35 

0.138 
0.136 
0.134 
§ 0.132 
S 0.130 

i 

0.128 
0.126 
0.124 

0.0 0.2 0.4 0.6 0.8 1.0 

a 

Figure 7. (Color online) Minimum average purity for n = 7 and an increasing number 
k of randomly selected bipartitions. Every point corresponds to an average over different 
extractions of k bipartitions. The bars correspond to 68% of the distribution (that can be very 
asymmetric due to the closeness of the lower bound), a = k/K. The average fraction of 
frustrated bipartitions is 1 — a c ~ 34%. 
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